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We consider general fermi-phobic scenarios in which excess events in diphoton or WW/ZZ res- 
onances may be seen at LHC. These Higgs like signals do not necessarily suggest that the new 
resonance is a particle with Yukawa couplings nor do we know that it is responsible for electroweak 
symmetry breaking. We can, however, extract two facts from it, this particle couples to pairs of 
SU(2) and U(l) gauge bosons and it must be a scalar, pseudoscalar, or tensor. We consider the 
signals of general operators up to effective dimension 5 in which a new scalar, psuedo-scalar, or 
O |. tensor particle may couple to pairs of standard model gauge bosons. This particle may or may not 

be charged under the standard model gauge groups, and may be produced via gluon fusion or EW 
vector boson fusion. 
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I. INTRODUCTION 



The Standard Model (SM) has been amazingly successful in describing fundamental interactions however, it remains 
incomplete. The Higgs boson is most likely the particle responsible for electroweak symmetry breaking, though it has 
yet to be discovered. The Higgs mass is not a known quantity however, its production cross sections and branching 
, fractions to SM particles are well known 1] 2] 3]. At the LHC both the production and the decay of the Higgs boson 
• are dominated by processes in which the Higgs has effective couplings to pairs of SM gauge bosons. While Higgs 
production is dominated by gluon fusion [U, it does not couple directly to gluons but has a large effective coupling 
them through a loop of heavy quarks. At low masses (m# < 130GeV), the Higgs decays predominately to bb however, 
t— I ' large QCD backgrounds at LHC make this channel difficult to see. Thus, the Higgs decays to pairs of electroweak 
gauge bosons are crucial for Higgs searches. In fact one of the most important channels for Higgs discovery is h — > 77. 
. i-H ' Again, the Higgs does not couple directly to photons but it has an effective coupling brought about by loops of charged 
W bosons and heavy third generation quarks. In fact, decays to weak boson pairs become dominate once we reach 
the kinematic threshold since the leading term in the partial width of Higgs to diboson is cubic in Higgs mass versus 
linear in Higgs mass for partial width to fermion pairs. 

LEP, Tevatron, and now ATLAS and CMS have been chipping away at Higgs mass parameter space leaving only 
a small window for the SM Higgs to hide [5(tS|- Recent search results from both ATLAS and CMS have reported a 
small excess in the WW channel and a slightly larger excess in the ZZ and 77 channels around 125GeV [7[|8(. While 
this excess may be the first hint at the SM Higgs, it is not yet a discovery. 

However, seeing a mass resonance in a diboson channel does not necessarily mean that the particle being produced 
is a Higgs. One would still be lacking evidence that the particle has the standard Yukawa couplings, or that it is 
responsible for electro- weak symmetry breaking. In fact, if a mass resonance is seen in and EW boson channel the 
new state may not have effective couplings to gluons as the Higgs does, unlike the Higgs whose production proceeds 
dominantly through gluon fusion, the new state may be dominantly produced through electro-weak couplings. A 
Higgs-like diboson signal may be the consequence of a variety of effective operators where a new state couples to pairs 
of standard model bosons. Such a particle may a scalar, pseudoscalar, or tensor; it may or may not carry standard 
model quantum numbers; and it need only have large enough effective couplings to SU(2) x C/(l) gauge bosons and 
may or may not couple to gluons. Various scenarios have been studied exploring some of these possibilities [9j floj [llj] . 

In this paper we present various Pseudo-Higgs scenarios in which a particle with effective couplings to SU(2) x U (1) 
gauge bosons may mimic a Higgs. We attempt to study general operators for all fermi-phobic scenarios up to effective 
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dimension 5. We explore constraints on the available parameter space and explore which regions of this space might 
be responsible for the observed di-photon excess at LHC. We pay particular attention to the possibility that new 
states may be produced dominantly by weak-boson fusion. The paper proceeds as follows, in Section 2 we consider 
the lowest dimension effective operators in which a scalar field may couple to pairs of SU(2) x U(l) gauge bosons- 
both the 'bosonic Higgs' scenario and a mimic scenario in which a scalar field couples pairs of gauge bosons through 
a modified Higgs kinetic term. In Section 4 we consider couplings of a SM singlet to gauge bosons through effective 
dimension 5 operators. We consider both pure scalar and pure pseudoscalar fields which are singlets under all standard 
model gauge groups. In section 5 we consider effective operators coupling a spin 2 field to pairs of SM gauge bosons. 




FIG. 1: Production of a SM singlet state through weak boson fusion 



II. REVIEW OF BOSONIC HIGGS 

Here we review the 'Bosonic Higgs', a non-standard Higgs scenario in which a particle is produced and decays to 
pairs of electroweak bosons thus mimicking a standard Higgs signal. In the model under consideration a field exists 
which has the quantum numbers of the standard model Higgs. The field is an SU(2) doublet which gets a vacuum 
expectation value breaking the SU{2) x 17(1) gauge group. However, this field does not have any Yukawa couplings. 
At tree level the Higgs-like particle couples only to W and Z bosons and the coupling occurs through the field's kinetic 
term. 

DH(DH)' = ((d-gA)H) 2 

We see that by expanding the covariant derivative one may take the piece proportional to A 2 and insert one power 
of the Higgs vev getting a term in the Lagrangian 

L ~ vHAA 1 

This is a term of effective mass dimension 3 coupling the Higgs-like field to the W and Z gauge bosons. Since there 
are no Yukawa couplings to couple it to heavy quarks this field has no effective couplings to gluons. The field does 
have an effective coupling to photons which it gets through loops of charged W bosons. 

Collider consequences of such a scenario have been considered by [13] [13] and recently by Q . 

In this case Higgs production proceeds dominantly through weak-vector boson fusion. For a Higgs with masses 
under 160 GeV, decay is dominated by the process h —> WW* /ZZ* . However the Higgs may decay to pairs of 
photons through a loop of charged W bosons. The Higgs decays to Z-f through a similar process. The total branching 
fraction of the Higgs to photons in this scenario is much larger than in a standard Higgs scenario. One reason for 
this is due to the lack of Yukawa couplings, this ensures that the branching fraction is not eaten up by decays to 
lepton pairs and to gluons. Another reason is because in the standard scenario, the Higgs couples to photons through 
loops of charged W's but also heavy quarks, which partially cancel the W loops. Without Yukawa couplings these 
cancelations do not take place and the partial width to photons is increased (see for example |l4(). Thus, though 
the total Higgs production cross section in such a scenario is smaller than through a gluon fusion process, the total 
production cross section of di-photons in this scenario may be appreciable, effectively mimicking a standard Higgs 
signal in the di-photon channel. 

A. New Effective Operators 

We now consider a scenario which has similar physics to that of the Bosonic Higgs case but with a much different 
particle. One may consider a particle S, which is a singlet under all gauge groups. One may write an effective coupling 
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FIG. 2: Plot of scalar mass vs. branching fraction to di-photons in the Bosonic Higgs scenario 

of this particle to the Higgs kinetic term which is not forbidden by any symmetry The operator is dimension 5 and 
suppressed by a scale A 

Here k is an arbitrary coupling. Now expanding the covariant derivative, extracting the A 2 piece we get 

^DH(DH)' = - gA)H) 2 -+ ^-AA^HH^ 

By inserting both Higgs vevs we get a coupling of the scalar field S to pairs of electro-weak gauge bosons. This 
effective operator is 

A 

has the same mass dimensions as the Bosonic Higgs operator but is multiplied by a factor kv/A. We see that a 
dimension 5 operator has become an operator of effective mass dimension 3. Also, in the Bosonic Higgs scenario, the 
coupling of the Higgs held to pairs of W and Z bosons was completely determined by the Higgs vev and the gauge 
coupling. In this scenario, however, the coupling to the gauge bosons may be varied by shifting the effective cutoff of 
the theory fc/A. As in the Bosonic Higgs scenario, a coupling to pairs of photons and to Zj will be generated at one 
loop through couplings to charged W bosons. When the scalar state decays, its branching fraction will be dominated 
by decays to pairs of W and Z bosons however, decays to di-photons will make up a few percent of its total branching 
fraction. Figure 2 shows a plot of scalar branching fraction into photons vs. the scalar mass for a particular choice 
of effective cutoff. The branching fraction of this state to photons is nearly 100 times that of a standard Higgs of the 
same mass. 

Like the Bosonic Higgs such a scalar state has no effective coupling to gluons. At LHC, such a state will be 
produced through the process involving electro-weak vector bosons, one such process is weak-boson fusion as seen 
in figure 1. Though this state would only have electro- weak production cross section, its relatively large branching 
fraction into photons means that it would have an observable signal. The entire process of production and decay would 
be qq — > Sqq — > jjqq. We may check to see if such a theory is consistent with LHC's observed di-photon excess at 125 
GeV. We simulated the production of di-photon events at LHC for 7 TeV with the process qq — > Sqq — > ~f~fqq. Events 
were simulated using MADGRAPH 5 [l5[ and showered with PYTHIA [l6j]. It was found that the acceptance for 
these events was similar to the acceptance of Higgs events decaying to di-photons and produced through the process 
gg -> h. 

Figure 3 a shows the overall production times branching fraction for di-photon events in our scenario vs. the variable 
effective cut-off scale for a scalar of mass 125 GeV. One sees that for sufficiently high cut-off the effective coupling of 
the scalar to SM gauge bosons is suppressed and the total production cross section falls off. The overall strength of 
the signal will depend on the scalar mass and on the effective cut-off of the theory. In figure 3b we present the region 
in the scalar mass plane vs effective cut-off plane and find it is consistent with the di-photon excess at LHC at 125 
GeV. 
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FIG. 3: Left- Contour Plot of total di-photon production cross section in the scalar vs. effective cut-off mass plane. Right- 
allowed region in effective cuttoff parameter space consistant with the LHC di-photon excess for a 125 GeV mass resonance 



III. DIMENSION 5 OPERATORS 



Alternately, one may consider operators of effective dimension 5 through which a particle that is a singlet under 
all standard model gauge groups may couple to many or all of the SM gauge boson. In this case the particle may be 
either a scalar or pseudoscalar. The Lorentz and gauge invariant operators one may write are 

L = hSF^F^/A 
L = hSpF^F^/A 

Where S is a scalar field and S p is a pseudoscalar. One sees that the scalar couples to the square of the field 
strength tensor, while the pseudo-scalar, which is parity odd, couples to the dual tensor. Unlike the Bosonic-Higgs 
like scenarios here the scalar or pseudo-scalar field may have independent couplings to each SM gauge group and thus 
there is an independent parameter for each gauge group. 

One may consider the parameterization 

K = k/A 

where we have absorbed the couplings into the inverse of the cutoff scales. As k* falls the effective scale of the cutoff 
increases. These three parameters, k* along with the mass of the new particle determine all of the new observable 
physics of the scenario. The coupling coefficients of the scalar or pseudoscalar particle to pairs of SM gauge bosons 
are given by: 



gsww — 


2k 2 


sw 2 A 


gszz = 


g\ ^sw 2 
4e 2 cw 2 


ffs 77 = 


g\ h + k 2 


4e 2 A 


gszj = 


.9152^ k 2 
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9Sgg = 


h 

A 



W )/A 



Many of these couplings may be varied independently. First we see that the SU(3) coupling is independent from 
all others. The new particles may or may not couple to gluons. The coupling to W's depends only on the SU(2) 
parameter k 2 . Thus, the coupling to W's may be turned off or varied while the coupling to photons remain the same 
by compensating with a change in k\. The couplings to photons and the Z boson may not be independently varied. 

Different values of the parameters allow for different modes of scalar or pseudo-scalar production to be dominant. 
If the couplings to gluons are at all appreciable then gluon fusion becomes the most important production mechanism 
for the new particle-it is produced through the process gg — > S. However, if the coupling to gluons is not present or 
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FIG. 5: For Left- pseudo-scalar and Right- scalar models with dimension 5 couplings to gauge bosons: regions in effective 
5(7(2) x (7(1) coupling space which fit the LHC di-photon excess for a 125 GeV resonance for varying values of effective 5(7(3) 
coupling. 



very tiny then the production process is dominated by the weak boson processes, qq —> Sqq. One finds in general that 
the production cross sections of scalar and pseudo-scalar particles are quite similar. 

In addition the branching fractions depend heavily on the choice of parameters. In general, the branching fractions 
into massive final states, WW* and ZZ* are quite small compared to that of massless final states 77 and gg. The 
branching fraction into di-photon may be large over a wide range of scalar or pseudo-scalar masses if the effective 
couplings kl, k* are big enough. If the coupling to SU(3) is zero or negligible then photons dominate the branching 
fractions. As the coupling to gluons is turned on then the branching fraction to gluons quickly becomes large. We 
have given analytic formulae for both scalar and pseudo-scalar branching fractions in the appendix. 

To illustrate the behavior of branching fractions, figures 4 a and b show contours where the photon branching 
fraction is 10 percent in the space of effective couplings fejf, k\. The contours are shown for increasing value of 
effective coupling k\. One sees that at small values of k*, the branching fraction to photon pairs will remain above 10 
percent even for small values of effective coupling k*, but as fcg-the effective coupling to gluons gets larger then 
the effective coupling to the weak bosons k* , k% must get much larger to keep the di-photon branching fraction above 
ten percent. 

Such a model with a scalar or pseudo-scalar may have Higgs-like signals. Throughout much of parameter space such 
a scalar or pseudoscalar may produce many di-photon events. Singlet production may proceed either through gluon 
fusion, or if the couplings to gluons are absent, through weak boson fusion. The entire processes are qq — > Sqq — ¥ jjqq, 
and gg — > S — > 77. In figure 5 we present the region of parameter space where a 125 GeV scalar or pseudoscalar 
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may be consistent with LHC's observed excess. Events were simulated using MADGRAPH 5 [15[ and showered 
with PYTHIA [16| while the branching fractions were calculated analytically. We have shown the allowed regions 
of parameter space both where the coupling to gluons is zero and where it is finite. One sees that for appreciable 
couplings to gluons, the overall production cross section of the singlet field is quite high, therefore in order to be 
consistent with experiment the di-photon branching fraction must be small, forcing one into a small corner of fcjf, 
parameter space. Note there is a large region of viable parameter space where the couplings to gluons is zero and all 
processes are electro- weak. 

Alternatively, one may look to the current Higgs searches to place bounds on the parameter-space of new scalar 
and pseudo-scalars, as these states may exist with mass values that have been ruled out for standard Higgs scenarios. 
For example one may translate the Higgs exclusions in the di-photon channel to upper limits on the production 
cross-section of di-photons for new resonances. Figure 6 for example gives upper bounds in k*, k% parameter space 
for a scalar particle with dimension 5 coupling to photons for a range of masses. 



IV. OPERATORS WITH TENSORS 



Another possibility for a particle which fake Higgs signals by decaying into gauge bosons is a spin 2 state. Such a 
particle should be a singlet under the Standard Model gauge groups, and of course must not get a vacuum expectation 
value. If we are considering operators up to dimension 5 which couple a spin 2 particle to gauge bosons we have two 
possibilities 



£ = h(r^ - 4T^)FPF„ p + ^(T^ - 2T^)m 2 A A^A u 

The first term is an operator of effective dimension 5 similar to those we have already considered. The second term 
is equivalent to the effective dimension 3 operators we have considered which may be generated by coupling the spin 
2 field to a Higgs kinetic term. These couplings are quite similar to those of a KK gravition to pairs of gauge bosons, 



see for example 17]. In the case of the KK graviton both sets of operators must occur together, but we may consider 
them separately. In addition, in the case of KK gravitions the coupling to pairs of different standard model gauge 
bosons must be universal, while in out more general case we are not under such a constraint. 

The first operator allows a spin two state to couple to all pairs of gauge bosons, gg, WW/ZZ, 77 and Zj. However, 
the second operator only allows the spin two state to couple to W and Z bosons, the coupling of the spin two state to 
photons must then proceed through a loop of W bosons similar to the bosonic Higgs. Because the effective coupling 
to photons in that case is of higher dimension, we will here only analyze in depth the phenomenological consequences 
of the first operator. 
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FIG. 7: For models of a tensor pseudo-higgs with coupling to SM gauge bosons: Left) Regions in effective SU(2) x U(l) 
coupling space which fit the LHC di-photon excess for a 125 GeV resonance with varying values of effective 5*17(3) coupling. 
Right) Contour plot of branching fraction into di-photons in effective SU(2) x ?7(1) coupling space for effective SU(3) coupling 
k* 3 = 1(T 5 

We see that the coefficients of the couplings to pairs of gauge bosons is the same as those given in section 3. The 
couplings to gluons and W pairs may be varied separately. Once again if there is no couplings to gluons then the 
production of the spin two state proceeds through weak boson fusion. In this case, though one has only weak-scale 
production cross sections, there is a large branching fraction into photons, therefore one may expect the spin two 
state to produce a significant number of diphoton events with a total cross section similar to that of the Higgs. 

As the coupling to SU(3) is increased, gluon fusion becomes the dominant production mechanism of the spin two 
state and the production cross section increases dramatically . However, once the tensor can decay into pairs of 
massless gluons they quickly begin to eat up the particle's branching fraction. In the appendix we have included some 
of the details of the calculation of the spin two state's branching fraction. Figure 6b gives a contour plot of the spin 
two state's branching fraction into photon pairs in kl ,k% space - the space of effective couplings to U(l) and SU(2) 
gauge groups for fixed value of SU(3)coupling. 

Such a state may fake a Higgs-like diphton signal. Figure 6a shows the region of parameter space in which a 125 
GeV spin two object may give an signal in the di-photon channel which is between .8 and 2 times that of a SM Higgs. 
The separate shaded regions show the allowed region for different values of effective coupling to SU(3). One sees 
that there is allowed parameter space where this coupling is and the production and decay of the spin two particle 
follows only from electroweak processes. One also sees that as the effective coupling to SU(3) is increased, the overall 
production cross-section of the spin two state increases greatly. In this case, in order to not overproduce di-photon 
events, the overall branching fraction into photons must be small, meaning that the allowed parameter space is in 
regions with small effective couplings to SU{2) x U(l) gauge groups. 

V. CONCLUSIONS 

We have shown that a variety of effective operators may produce Higgs-like di-boson signals at LHC. Scalar, pseudo- 
scalar and tensor particles which are uncharged under SM gauge groups may couple to SM gauge bosons. The new 
effective operators allow for states to be produced either through gluon fusion or processes such as weak boson fusion. 
One may use current Higgs searches to place bounds on the parameter space of these models, in addition, these new 
states may have a large enough production cross section and branching fraction into di-photons to account for the 
125 GeV excess in the di-photon channel at LHC. 

Further areas of study may involve building completions of these models, such a task may or may not be difficult to 
achieve depending on the scenario. In the case of dimension 5 operators which couple a scalar or pseudoscalar state to 
the SM gauge bosons, completions may simple. The scalar may have couplings to TeV scale fermions charged under 
the SM gauge groups. Spin two state may be new light composite states. Any completion would likely involve new 
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particles between 1 and 10 TeV which interact with the SM gauge groups and should have experimental consequences 
for searches conducted at 14 TeV. 
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B. Branching Fraction Calculations 

1. Standard Model Singlet 

The 3-point interaction from ([lj that contributes to the decay of the scalar has the form 

C = (dpAvd^A" ~ d»A v d v A») 

This vertex has the Feynman rule —igsA 1 A 2 {pi ■ "Pig^u — Pi^P2v), where pi and p2 are the momentum of the gauge 
bosons. 

The decay width into massless states is straightforward to calculate. We find 

r(5^77) = 



47T 

r ( ,v-. „<,) = M^rn 3 



z 



The decay into Ws and Zs is a little more involved. At least one of the massive gauge bosons will be off shell and 
decay into fermion pairs ala [llj]. The matrix elements for these two cases are 

where p\ and P2 are the momentum of the fermions, P3 the momentum of the final state gauge boson, and q = p± +P2 
is the momentum of the off shell gauge boson. In terms of invariant mass, the decay width to massive gauge boson 
pairs becomes 



nS^Zff) = l 9szz 2 T 3 Z [ dml 2 dm 2 23 \Mf 
lob-K m%mz J 

T(S^Wff) = gg|gg T 3 W fdm\ 2 dml,\M\ 
1287r z vn%raw J 
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M < 2 = 71 2N2,r2 2 ^ m \ ^ [( m l - m i2 -m| 3 )(m| 3 - m 2 A ) +m 2 A m\ 2 ] 

\Q —m A ) -\- 1 " 



A" l A 



in 



2. Pseudoscalar 



The pseudoscalar fake Higgs can decay into pairs of gluons, Z's, W's, photons and to Z plus photon. The piece of 
the Lagrangian which contributes to these decays has the form 



Se^d^dpA^ 
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The Feynman rules for these vertices are given by —igsA 1 A 2 € ^' uf " T PipP2p where the prefactors are again given in Section 
3. The decay width to massless pairs and to Z plus photon are straightforward to calculate. They are 



47T 



r(s->z 7 ) = ^ 

7T V 



For the WW and ZZ channel, at least one of the massive gauge bosons will decay to a pair of fermions. In terms of 
momentum, the square of the matrix elements are 

"W W " l W ' ' W W 

\M\% Z . = 12^ (g2 _Jr +r|m| [2, 2 (, • - - 4tfV • ■ PS)] 

We do the 3 body phase space integral in terms of invariant masses, mi 2 and 11123- For both cases, the square of the 
matrix element is just a polynomial in m| 3 and so this integral can be done analytically. Because of the limits of 
integration on the to 23 integral, we do the integral numerically. In the rest frame of the decaying particle, the 
decay widths into Wff and Zff are 

r.2 



m? 



M = 



'12 



W 2 -^) 2 +rK 

and is either the W or Z mass. 



1/ 2 2 2 \2 / 2 2 2 w 2 2\ 22 

-(m s -m A - m 12 ) - (ra s - m 12 - m 23 )(m 23 - m A ) - m A m 12 



3. Tensor 



The completeness condition for the polarization tensors for T^„ is 

s 

The spin two pseudo-higgs can decay into pairs of gluons, photons, Zj, WW*, and ZZ*. The massless case and the 
Z-f case can be calculated analytically. We find 

T(S -> 77 ) = 
r(5^. 9ff ) = 2 -^m 3 s 

r,2 



/32m| 


22 m| 


8 


V 3 mi 


3 m* 





fci • fc 2 = m s \ m z 7 + m 2 s +m 
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One of the massive gauge bosons will decay into fermion pairs for the massive diboson case. The matrix elements 
squared in this case are 

m w (q 2 -m z w y + T z w m z w \ rn z w J 

.as ~\~ ^mn,as) 

= *% v-M<i<4 im - f " ■ l8 " + * * (-*" + If) 

where ki and &2 are the four momenta of the final state fermions, &3 is the four momentum of the final state gauge 
boson and q is gauge boson propagator momentum. The index contraction for these matrix elements was done using 
FORM [? ]. The three body phase space integral, done in terms of invariant mass, was computed numerically. 



S. Dittmaier et al. [LHC Higgs Cross Section Working Group Collaboration], arXiv:1101.0593 [hep-ph]. 

S. Dittmaier, S. Dittmaier, C. Mariotti, G. Passarino, R. Tanaka, S. Alekhin, J. Alwall and E. A. Bagnaschi et at, 

larXiv:1201. 3084 [hep-ph]. 

S. Dawson, |hep-ph/ 9901280 

J. Baglio and A. Djouadi, JHEP 1103, 055 (2011) |arXiv:1012. 05301 [hep-ph]]. 

[TEVNPH (Tevatron New Phenomina and Higgs Working Group) and CDF and DO Collaborations], larXiv: 1203.37741 
[hep-ex] . larXiv:1202.T488l [hep-ex] . 

R. Barate et al. [LEP Working Group for Higgs boson searches and ALEPH and DELPHI and L3 and OPAL Collaborations] , 

Phys. Lett. B 565, 61 (2003) hep-ex/0306033]. 

[ATLAS Collaboration], larXiv: 1202.14081 [hep-ex], 

S. Chatrchyan et al. [CMS Collaboration], 

E. Gabrielli, B. Mele and M. Raidal, larXiv:1202. 17961 [hep-ph], 

J. Ellis and D. S. Hwang, arXiv:1202.6660 [hep-ph]. 

I. Low, J. Lykken and G. Shaughnessy, Phys. Rev. D 84, 035027 (2011) arXiv:1105.4"587l [hep-ph]]. 
A. G. Akeroyd, Phys. Lett. B 368, 89 (1996) [hep-ph/95ll347] 



B. Abbott et al. [DO Collaboration], Phys. Rev. Lett. 82, 2244 (1999) hep-ex/9811029 

D . Rainwater, |hep-ph/0702124| [HEP-PH] . 

J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer and T. Stelzer, JHEP 1106, 128 (20 11) |arXiv:1106.0522l [hep-ph]]. 

T. Sjostrand, S. Mrenna and P. Z. Skands, JHEP 0605, 026 (2006) hep-ph/0603 175]. 

T. Han, J. D. Lykken and R. -J. Zhang, Phys. Rev. D 59, 105006 (1999) |hep-ph/9811350] , 

T. G. Rizzo, Phys. Rev. D22 , 722 (19 80). 

J.A.M Vermaseren, |math-ph/0010025 



